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1 Introduction 


The equations that describe the propagation in a Raman active medium 
when Stokes Eg, anti-Stokes Ea and pump Ep exists and when there is no 
frequency mismatch can be written B I, I, I, I, I in the form: 


dEp 

dC 

dEa 

dC 


= (3aQ*Ea-QEg, 


— ~l3aQEp, 


dEg 

dC, 

dr 


= Q*Ep, 


+ QQ — E*Ep + l3aE*Ea, 


( 1 ) 


where Q is the normalized effective polarization of the medium, ( = z/L 
and r = t — z/v are dimensionless space and retarded time coordinates, 
respectively, v is the wave group velocity. T 2 is the natural damping time 
of the material excitation and g = I/T 2 . By jSa we denote the coupling 
coefficient which determines the number of anti-Stokes photons relative to 
number of Stokes photons and its magnitude depends on the matrix element 
that describes the dipole transition 0 0- In this paper we consider f3a = 
1 B) but some results are valid for (3a ^ 1- For (3a = 0 i.e. when the 
anti-Stokes wave Ea is neglected, these equations are the so called transient 
stimulated Raman scattering (SRS) equations, which possess ZS-AKNS B Q 
representation when g = Q [T^. The SRS soliton solutions theoretically 
discovered by Chu and Scott have been experimentally observed first in 
pd]] - The SRS solitons (regarded as transient solitons |^, with a vr 

phase jump at the Stokes frequency) have been extensively studied [^. In 


later experiments by Duncan et al a careful comparison between theory 
and experiment showed good agreement. Shortly after this work |]T^, Hilfer 
and Menyuk [|^ carried out simulations which indicate that in the highly 
depleted regime the solutions of transient SRS equations always tend toward 
a self-similar solution. This result has been recovered by Menyuk et al [p!B|, 1^ 


applying the inverse scattering transform method (ISM) to the transient SRS 
equations, see |T2|. Experiment to observe this solution has been proposed 


in 1^. The Kaup’s theory |]^ also indicates that the dissipation, which 
appears for hnite T 2 plays a crucial role in soliton formation. The similarity 
solutions and other group invariant solutions of the SRS equations in the 
presence of dissipation are studied in ||2^ . Claude and Leon reformulated 
the transient SRS equations as an equivalent c)-problem and thus were able 
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to treat the inhomogeneous broadening with initial conditions of Driihl et 
al [M experiments. As a consequence they showed the that Raman spike 


observed in the experiment is not a soliton. 

For general system (|ID with anti-Stokes wave, phase mismatch, Pa ^ 1 and 
dissipation, transient tt solitons have been investigated by M. Scalora et al 0 
using numerical methods. They predict the formation of solitonlike pulses at 
the anti-Stokes frequency. Another type threshold bright 2 n solitons, which 
have Lorentzian form have been theoretically obtained in . 

In sect. II we introduce new variables S 3 , S± (@) which are quadratic in 
terms of Ep, Eg, Ea- Then the system (D for S 3 , S± derived from (|l|) with 


/da = 1, ^ = 0 allows ZS-AKNS representation similar to the one used by 
Chu and Scott for another physical quantities: difference of normalized 
Stokes-anti-Stokes local intensities and for normalized (complex) local Rabi 
frequency. We also introduce “nonlinear time” and renormalized dimension¬ 
less variables different from yhe ones used in [^. Then we solve the inverse 
scattering problem (ISP) for the system (PT|) with the “nonlinear time” r' 
restricted to the finite interval 0 < r' < 1; i.e. we derive the corresponding 
Gel’fand-Levitan-Marchenko (GLM) equation. 

In sect. Ill we obtain new periodic, soliton and self-similarity solutions of 
Stokes-anti-Stokes SRS equations without dissipation are obtained for both 
formulations (|^) and (|^). In addition the transient solitons and the bright 
solitons of Kaplan et al 


22 are discussed. 


In the last Sect. IV we propose an extension of the Stokes-anti-Stokes 
SRS equations for N Stokes and N anti-Stokes waves and conjecture that it 
is also integrable by means of the 1ST method. 


2 ZS—AKNS Representation and GLM Equa¬ 

tion 

2.1 ZS—AKNS Representation 

Let us introduce the following variables 

S3 = i(|s,l"-|s.t), s^='-{e:e, + e;e^), s+ = s:. (2) 
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In terms of the new quadratic variables (j^) the initial system with g = j3a = 
1 is rewritten as 

^*^3 dS+ dQ 

QQ *S^+ “1“ iQB—^ 2^o_i_. (3) 

Then the eq. (|^) can be written down as the compatibility condition 

drU - diV + [U, V] = 0, (4) 


of the following linear systems: 


L{X)F{C,t,X) ^ —-U{C,t,X)F{C,t,X) = 0, (5) 

M(A)F(C,t,A) = ^-l/(C,T,A)F(C,r,A) = F(C,r,A)C(A), (6) 


where 


U{C,r,X) 

qiC^r) = 


+ ^qiC,r), v{c,T,x) = 


S{C,t) = 




( 7 ) 

( 8 ) 


and C'(A) will be hxed up below. 

From physical point of view BH the initial value problem associated 
to the system (|I|) is the following 


Q(C,0) = 0, Ep{0,T) = Epoir), (9) 

E,(0, r) = Eso{t), Ea{0, r) = Eaoir) 


and the problem is to determine the output quantities Ep{L,T), Es{L,t), 
Ea{L,T), where L is the total length of the beam path in the Raman cell. 
Analogically the initial value problem for system (|^) is 

Q(C, 0) = 0, ^ 3 ( 0 , r) = ^ 3 o(r), ^+(0, r) = A+o(r). (10) 

We follow the main idea of [^, namely that as a Lax operator one 
should consider the operator M(A) in (P) and solve the inverse scattering 
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problem for it. Then we will use the second operator L{X) in and deter¬ 
mine the ^-dependence of the corresponding scattering data. However there 
will be substancial diferences in the details. 

First of all we will approach the inverse scattering problem for the M(A) 
operator directly rather than via its gauge equivalence to a ZS-AKNS type 
system. Indeed, this equivalence is realized with A) evaluated at A = 

0. However, in our case the other linear problem has a pole singularity at 
A = 0. Due to fact makes one can not evaluate the C^dependence of the gauge 
function, which makes impossible the comparizon with the results in |^, ^ . 

It is well known how to solve the ISP for the system considered on the 
whole r-line —oo<t<oo and with boundary conditions of ferromagnetic 
type, i.e. limT-^±oo r) = see 0. We will make use of these ideas 
adopting them to our case. First we have to take into account that the 
eigenvalues of our S{(, r) differ from ±1 and are generically r-dependent. In 
order to calculate them it is enough to know, that tr ^(C, r) = 0 and 


-detA(C,r) = Sl + 2S+S_ 

= i ( l^y - + \\e;e, + e,e:\^ = k*(t) 

Using the evolution equations (|l|) we check that 

dK 


( 11 ) 


dC 


= 0 . 


( 12 ) 


Besides from ([Tl|) we conclude that K{t) is real-valued function. In order 
to proceed further we require in addition that is monotonic function 

of r. Then we can introduce a new “nonlinear time” t' by 


dr' = K‘^{T)dT 


(13) 


and the following dimensionless variables: 


r' = / K\r")dr"lT^, 
Jo 


T = 

OO 


pi ^ 

P KirY 


E' = 




poo 

/ KYT)dT, C' = CToo (14) 

Jo 

Q 


EL = 


E. 


Q’ = 


K^tY “ K{tY 

The primed variables introduced above satisfy the same system (13) of NLEE 
provided g' = gToo/K^; in what follows below we put ^ = 0. Note, that the 
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transformation {Ep, Eg, Ea,Q} {E^, E'^, E'^,Q'; K{t)} is one-to-one and 
invertible. In order to obtain the evolution of {Ep, Eg, Ea, Q}, one must first 
determine the evolution of {E'^, E'^, E'^, Q'} and then use the given function 
K{t) to return to the original variable set. 

The nonlinear time t' is introduced in analogy to the one in ; the dif¬ 
ference is that now K^{t) can not be interpreted as the total energy density: 


S{T) = \Ep{C,T)\^+\Eg{C,T)\^+\Ea{C,T)[^ 


(15) 


which is constant at every point r as function of Since all physical solutions 
possess hnite energy we conclude, that each of the terms in ([l5|) must be 
integrable functions of r. In particular, each of these functions must vanish 
for r —cx). As a consequence of this fact we conclude that K^{t) must have 
the same properties. Therefore for this class of solutions we have T^o < oo 
and in terms of r' we get the system: 

M'(A)F(C, t', A) = - AS'(C, t')F(C, t', A) = 0, S'((, r') = 

lie) 

where S\(, t') satishes tr S" = 0 and det S' = —1. As a result the eigenvalues 
of become equal to ±1, i.e. we can write down 

S'iC, r') = g{C, r')(T^g-^{C, t') (17) 


Note that the constancy of £{t) corresponds to pointwise conservation 
of the photon intensity. From Eq. ([^) it is easy to derive also the following 
important relation: 

+ ^*^3 = 0 . ( 18 ) 

from which for ^ = 0 we hnd that / IQpdC is an integral of motion if 
>S'3(0, r) — S' 3 (L, r) = 0. This will be fulhlled if Ea and Eg satisfy quasiperi- 
odic boundary conditions, i.e., if Ea^gi^r' = 0) = Ea^sij' = 1) with any 

<Pa,s- 

In the next subsection we will use only renormalized quantities and the 
“nonlinear time” t' and for the simplicity of the notations will drop all primes. 
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2.2 The GLM equation 


We briefly sketch the derivation of the GLM equation related to the left end 
r = 0 of the interval. Of course we have to introduce also slight modihcations 
in order to take into account the fact that S'(r = 0, C) = 5'o(C) 7^ o's- The op¬ 
erator M(A) on hnite interval generically possesses purely discrete spectrum 
with an inhnite number of simple discrete eigenvalues. As a consequence, the 
kernel of the GLM equation contains only a sum over the discrete spectrum. 
Skipping the details we write down the results. 

Let the dost solution of ([T6|) normalized to the left end r = 0 of the 
interval, be hxed up by: 


0o(gC,A) = lim0(r,C, A) = lim0o(C, A) = 11, (19) 

r —^0 r —^0 

= limgi(C)e"'">'“r(C,A), (20) 

r —^0 r —^0 

^o(C) = ^(C,^ = o), giiC) = g{C,r = i), ( 21 ) 


Then its behaviour at r ~ 1 determines the scattering matrix T{(, X) accord¬ 
ing to (^). It remains now to evaluate the “evolution” of T in (. In order to 
do this we have to calculate hrst G(A) in (|^) by taking the limit of ( |T^) for 
r —> 0 with the result G((C, A) = (V'^) 5 'o"^'^ 35 'o(C)- Then we take the limit of 
(P!9[ ) for r —A 1 which gives the following result for the evolution of T{(, A): 

^ J (r(C. A)9„-'<T3g„ - gr'<T39ir(C, A)) (22) 

As we already noted, in order to apply the ISM for the solution of our problem 
we will need hrst to calculate not only goiC), which is determined from the 
initial conditions, but also gi{C)- The situation is greatly simplihed if we 
impose a quasiperiodic boundary conditions on the helds Ea^s^p{(, r) in such 
a way, that S'o(C) = *S'i(C). Then we get go(() = gi(() and the r.h.side of (^) 
becomes proportional to the commutator [T, 5 'cr^(C)c’' 35 'o(C)]- The importance 
of this imposition can be seen from the fact, that it immeadiately provides 
us with the hierarchy of conservation laws. The generating function of this 
hierarchy is trT(A), which is now (^^independent. 

Let us also be given S'o(C) = = 0) and let it be diagonalizable in 

the form: 

SoiC) = go{C)(^3go^iC), ^o(C) = = o). (23) 
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We introduce the transformation operator which relates the dost solution 
0(r, (, A) to its asymptotic (, A) (|^): 

A) = 0o(C, A) + / r_(r,^;C)0o(^,C, A)ci^ (24) 

Jo 

We have to keep in mind also that all solutions and the scattering matrix of 
the system (p!6D are meromorphic functions of A. 

Then we obtain that r_(r, ?/;C) must satisfy the following GLM-type 
equation: 

^-ir,y;C) +So{C)K{T + yX)+ f T_{t, z;C)K\t + z'X) dz = 0, (25) 

Jo 

where the kernel K{t] () and its derivative K' = dKij] O/dr are given by: 

K{r-, 0 = 9„(C) ( . 0 ) 9„-‘(C). Mr; C) = - E (26) 


Here Xj are the discrete eigenvalues of M(A) and mj{Q is related to the norm 
of the corresponding dost solution of ([l6|) ; generically Xj may depend also on 


C- 

The corresponding potential of ([I^) is recovered from the solution r_(r, ?/; Q 
of (P^D through: 


S{T,C) = B.{r,C)So{OBz\T,C), 


i?_(r,C) = Il + r_(r,r,C)5o(C). (27) 


The complete solution of the problem requires also the calculation of the 
C-dependence of the scattering data, in our case rrijX) and Xj. 


3 Periodic, Soliton and Similarity Solutions 


In this section we generalize the results of [O, 20 to describe the similarity 


solutions-solitons, periodic and self-similar solutions for the system (^p. 

We also have not been able to resolve the fundamental problem, inherent 
to this type of NLEE. We have to solve the ISP for M(A) on a finite interval, 
and naturally the (^-dependence of the corresponding scattering data will 
depend on the boundary values of Q at both end of this interval. Indeed, 






the initial conditions allow us to calculate all the necessary quantities such 
as r = 0), S±{(, r = 0) at r = 0. In order to evaluate them at r = 1 
we have to solve completely the problem. 

On the other hand, the initial conditions of such physical system must 
determine uniquely its evolution. One way out of this problem is to im¬ 
pose certain boundary conditions on the operator M(A) - e.g., (quasi-) pe¬ 
riodic. They will relate the values r = 0), S±{(,t = 0) to r = 1), 
S±iC,r = l). 


3.1 Periodic (cnoidal) and solitary waves 

At hrst we will study the cnoidal wave similarity solutions which include 
solitons as a special limit. If we introduce the retarded coordinate ^ = (—r/a 
and the following new variables 

BAO = t), B.iO = r), 

= y(C) = e-“""Q(C,T), ( 28 ) 


where e, ep^s,a are arbitrary real parameters. We hnd that the reduced vari¬ 
ables Bp, Bg, Ba and Y satisfy: 


dBp 

dBa 

di 


PaBaY*e-^^^ - e^^^BgY, = e-^^^BpY*, (29) 

1 

-fd^e^^^BpY, -— + e-^^^BpBl + l3ge-^^^B;Ba = 0, 

a d^ ^ 


where ep = Eg + e, -|- e. The general solution to Eq. (|^) can be 

expressed in terms of elliptic integrals. When e = 0, the solutions may be 
explicitly written in terms of cnoidal functions. Indeed, from (^ and (|^) and 
using the transformed variable ^ — r/a we hnd the system: 




-iQ*S+ iQS-, 


dS+ 


-iQSs, 


dQ 


which has the following hrst integrals 


(30) 


-SJ + S+S_ = h, S+Q- + S.Q = I 2 , 


-S3 + i|<3P = /3. (31) 

a 2 
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Introducing the new real variables A^, Q and 0 by: 


S+ = e^^+A+, 


S- = e 


Q = 


= 0-0+. (32) 


we rewrite (^) and (^) as follows: 

dSs „ ~ ~ dA^ 


d^^ 

dQ 


a 


= —2QA_^_ sin 0, 

= 2^4+ sin 0 , 

1 


h = 2^3 + 


d^ 

dcj) 2 A^ 

d^ a Q 

I 2 = 2 QA^ cos 0 , 


= (55'3sin0, 
cos 0 , 


(33) 

ds = —S 3 + -Q^- (34) 
a 2 


Squaring the equation for S 3 and using the expression for fc = 1, 2, 3 we 
obtain 

42 ^ - Z 2 )iS 3 - Z 3 ), (35) 


a 


where the constants Zj are related to Ik by: 

Zi + Z2 + Z3 = a/3, Z1Z2 + Z2Z3 + Z3Z1 = —2/1 
Z 1 Z 2 Z 3 = a/|/4 — 2a/i/3. 


(36) 


The solutions for S 3 may be written explicitly in terms of Jacobian sn 
function. We have the following periodic (cnoidal) waves: 


1. For a > 0, Zi < 0 < Z 2 < Z 3 , 

53 = Z i + (Z 2 - Zi)snXe - ^o), k)], 

'3 — Zi 2 Z 2 — Zi 


(37) 


P = 


= 


a 


Z 3 — Z\ 

2. For a < 0, Zi < Z 2 < 0 < Z 3 , 

5'3 = -^3 + {Z 3 — ■Z^2)sn^[p(^ — ^ 0 )) k)], 
IZ3 — Zi 2 Z3 — Z2 


(38) 


P = 


= 


—a 


Z3 — Zi 


In the particular case when A; = 1 we find the corresponding solitary 
waves 
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3. For a > 0, Zi < 0 < Z 2 = Z 3 , 



(39) 


4. For a < 0, Zi = Z 2 < 0 < Z 3 , 


S3 — Z2 -h 





(40) 


Let us concentrate on the most important from the physical point of view 
soliton solutions, i.e. the third case with additional constraint Z\ = Z 2 = 
—aJa- The result of integration is 



(41) 


where ^0 is the arbitrary initial phase. We will return again to this solution 
in section 3.4. 

3.2 Self—similarity solutions 

In order to obtain the self-similarity solutions we introduce the reduced vari¬ 
ables 



where ^ = ^r. Then the equations (J^) becomes 
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(43) 













where €„ = + e, ea = €„ + e. For e = 0 


simplifies to: 


= ~B Y + -5 Y* 

di t‘ r° ' 


^ = -B Y* 


dBa fd 




= -jBpY, 


dY 

— = b,b: + (3b^b;, 


(44) 


and, while the general solution is singular, nonsingular solutions in terms of 
series can be obtained by a technique described in 10 ■ 

We prefer here to analyze the solutions of with another self-similarity 
variable ^ = 2y/2(T. If we choose 


S 3 = cos[;3(0], 'S+ = J v'2sm[/3(0] 


(45) 


we find that 


goes into: 

d^m , ^dm 




+ - 


e d^ 


+ sin[^(0] = 0, 


(46) 


Equation (^Gf) was first derived in another context for the transient stimulated 
Raman scattering by Elgin and O’ Hare This equation can be reduced to 


one of the standard forms of the Painleve (-Pm) equation When S> 1, 
we can use the asymptotic formula given by Novokshenov to obtain 


/9(0 = 7 ^ cos (.^ + ^ In ^ + -0 


^ 1/2 


16 


(47) 


where 

16 

=-ln[cos(/3o/2)], 

71 

21ri2 / 7 (y‘^\ tt 

ln[cos(/?o/2)] TargF f — j - (48) 


Here r(a;) is the Gamma function with a complex argument, arg[r(a;)] indi¬ 
cates its phase and Pq = = 0). Similar expressions can be obtained for 

and 5+ from (|45D and for Q from 
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3.3 Discussion 


To obtain the bright solitons and to compare onr resnlts with the ones of 
we slightly generalize eqnations (|I]) (see for example 


Kaplan et al 


0 ). The Raman qnantnm transition between the lower (gronnd) and upper 
(excited) level, i.e. two level atom is described by a 2 x 2 hermitian density 
matrix p and the generalized Bloch equations 


dr ^ « ’ 


5A 

= 2 Re(( 5 f 2 /j), 


— 


^ ( 
h v" 


+ a 


s,p-^s 


E E' 

p,a-^p-^c 


^ — Pn — P 22 i 
Q = -2*^126-*^°^+*^"°'. (49) 

CUo/c 


hn hrt kc 


Here Clji is the generalized local Rabi frequency p6 
and we have assumed that pn + P 22 = 1- 

The generalization of (^ we mentioned above is obtained by replacing 
the equation for Q in (|l|) by the system (^9|) . These equations, rewritten 
for the self-similarity variable ^ = ( — r/a coincide with equation (6) from 
Kaplan et al. [0. The direct comparison of and Clji shows that the 
physical interpretation of S'+ is the normalized local Rabi frequency. From 
these equations we obtain also that the quantity 


A^ + |Qr = C(C), 

is conserved in r and that the following equations 


(60) 


— Ht), 
d 


\V„ 




d 


i = p,s, a, (51) 
(52) 

hold. By Nq we have denoted the density of Raman particles. Let us use the 


_(5a<ha-5.<hs)-7rAo—A = 0, 
oC or 


Vn 


<h = \E,, 


following ansatz [22 


— |ap| 'f’s, 4>s — [Osj d’a — jOal *^2, 


la P = - 


laJ^ = 


w 


s,pi 


52hF ’ 


\an? = 


\w. 


p,a \ 


27r 


Ws,p = 


-a 


s,p\ 


IUJgiOp 
^s^p 
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27r 


'^p,a 


-a. 


p,ai 


52 IF ’ 

iLJpLtJa 
ripTla ’ 


(53) 























where 0 


m 


1 


lm*^s )(h m2-6p) 


(^ml ^p) 


(<^ml + 


(*^lm*^a) (^m.2*^p) ^ (^lm*^p) (^m2*^a) 


fl m (^ml ^a) (^ml “(“ ^p) 


(54) 

(55) 


Here dim, <^ 2 m are the dipole matrix elements between the Raman quantum 
levels and the m-th quantum level, m ^ 1,2 and np,ns,na are refractive 
indexces at frequencies ujp,uJs,oJa respectively. From (|51D, we hnd that 
if ujp S> uJa, 0 Js and oJa,ojp ujmi then Ug^p — ap^a and ujg^p — oJp^a- Using 
arguments analogous to the ones in we hnd that one can expect such 
physical systems to be described by the system (|) with ~ 1. 

) we obtain: 


Inserting 


into / (1^ , from 


73 


\w 


s,p I 


dgSp 


+ 


m 


p,a I 


daSp 


Let us consider the second integral (|5^) with conserved density 5a<Fa— 
and conserved hux- ttNoA. Then 


Here 


T o r .T. 5 27r<FsfVo 

J = 2 (5a$a ^6*^^) = . 

<±>0 

(56) 

AttNoW 

° {\WsJ^/6s - \Wp^a\‘^/6a)' 

(57) 


where ” indicates that the molecules (atoms) are initially at the equilib¬ 
rium and ”-l-” that the population difference is inversed. Let us also intro¬ 
duce: 


d>s = $o^(0, A = ±(1-R(0), Q{0 =-73^3(0, (58) 

where S{^) may have different forms (soliton, Lorentzian etc.) Finally 
from the normalization condition |app -|- -|- joap = 1 we have 

|2 


w = 


w 


s^p \ 




p,a I 




dadp^a 
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(59) 

















where 


5 


s,p — 




(60) 


Recently these bright solitons, in more general physical sitnation, cascade 
SRS 1^, 1^ have been nsed to predict generation of snbfemtosecond coherent 
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pulses in SRS experiment 


From the above analysis it is clear, that 


bright solitons are obtained in the case of hnite group velocity dispersion 
parameters Si (0 i 1^ and non-zero population difference A. 


3.4 The one soliton solution 


Here we will show, that the auxiliary linear problem for the vector NFS 
equation with some additional reduction is equivalent to the Stokes-anti- 
Stokes SRS equations without last equation for Q. This formal equivalence 
allows us to recover Ep, Eg, Ea from potential Q, already obtained by inverse 
scattering transform method of Sec. 2 with (|T^. 

Indeed, if we introduce: 




( -01 
^2 
\ V’3 


/ 0 

Qi 

(l2 \ 

( 

-Qh 

0 

0 

^2 

V -?2 

0 

0 ) 

Ws / 


(61) 


and require that: 


gi = -Q, q2 = Q*, V’l = Ep, = Eg, V’s = Ea. 


(62) 


without last equation for Q, which we may obtain from the previous sub¬ 
section. Using the well known one soliton solution of the vector nonlinear 
Schrodinger equation under the reduction (|6^ and the solution of the linear 
problem (|UT|) with potentials Q, Q* we obtain 


Q = 


\/2rjt 


cosh( 2 ;) ’ 


z = riC- 


T] Jo 


K‘^{T')dT' 


(63) 


where K{t') is real, the soliton’s eigenvalue is irj and 0 is constant real phase. 
The direct integration of Eq. (^) with reduction (|62D is given by 


Ep = V2 K{t) 


tanh(^) 

cosh(^) 


e*'^, Es = iF(r) tanh^(z), Ea = 


K{r) 

cosh^(^) 




(64) 
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These solutions are similar to transient SRS solitons obtained in [^] (see also 
and for it' = 1, 0 = 0 coincide with the ones in [|I|. If we now calculate 
the S 3 and S± nsing the above expressions for Ep, Eg, Ea we obtain precisely 
the soliton solntion (|4T|). 


4 A^—component generalizations 


In this section we show that the extended model with A^-Stokes and N- 
Stokes components is also integrable in the sense means of 1ST method. The 
considerations are formal from physical point of view. 

Let us consider the following eqnations 


” = Y.WEf - QEf), = Q-E„ 

N 


dC .=1 

= -PaQEp, 


i = 1, 




dC 

^ + 9Q = E(d‘’'£'p + i3.e;e<:>). 

i=l 


(65) 

with /3a = 1. The eqnations for the E^ can be written down as the same 
auxiliary linear problem, which solves the A^-component vector NTS eqnation 


Qi — —Qi 0.2 — Q* I ■ ■ ■ (I2N-1 — —Qi (l2N — Q*-, 

= ... <p 2 N = Ei'^\ = 



f A ^ 


( 

0 

Ql 

92 • 

• 92 V-1 

92V 

\ 




^2 




0 

0 . 

0 

0 




d 

^3 



-<l2 

0 

0 

0 

0 



"03 


'lp2N 




0 

0 . 

0 

0 



'02V 


V V’2V+1 ) 


\ 

~^2N 

0 

0 . 

0 

0 



V 02V+1- J 


with the spectral parameter A = 0. 

We introduce again the quadratic variables 


1 JV ■ N 

S3 ^ - \E.W s+^si = -J2(e:e, + e;e,), ( ee ) 
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and show that they satisfy the same equation (|^). Therefore the ZS-AKNS 
representation (|) can be used as above for the analysis of the system (| 6 ^) . 
The procedure of solving Eq. (^) is analogous to the considerations of 
section IV. Clear physical interpretation and solutions of Eq. (|65D will be 
given elsewhere. 


5 Conclusion 


In this paper a method of solving Stokes-anti Stokes SRS with = 1 is 
presented. New transient solitons (^), (|64|) are obtained. For bright solitons 
our results are in agreement with these of Kaplan et al [^. The traveling 
and self-similar solutions to Eq. (|l|) are discussed. 

We stress also that for the transient solitons with f3a = 1 the number of 
Stokes-anti-Stokes photons are close to each other. 

The ISM generically meets with difficulties due to the fact that generically 
no boundary conditions on the potentials ([T 6 |) are imposed. In the case of 
quasiperiodic boundary conditions these difficulties are overcome. 

Looking at the Lax operator (pT]) we recognize the system of equations 
@ as belonging to the Heisenberg ferromagnet hierarchy. One can try to 
apply the expansions over the “squared solutions” in the spirit of and 
then treat the case ^ 7 ^ 0 as a perturbation. 
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